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In this study, a prototype helium-changed twin-power-piston y-type Stirling engine has been built, and 
some of its geometrical and operational parameters have been investigated by a numerical model. Data 
taken from the prototype engine have been used to correct the values of some factors in the numerical 
model. The results include volume and temperature variations in the expansion and compression cham¬ 
bers, p-v diagrams, and the effects of regeneration effectiveness, the crank radius of the power piston, the 
initial pressure of working gas, and the rotation speed on engine’s power and efficiency. It has been found 
that regeneration effectiveness poses the most prominent effect on efficiency, while engine speed is most 
effective on the engine power within the range of engine speed investigated in this study. This study 
offers invaluable guides for the design and optimization of y-type Stirling engines with similar 
construction. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

The Stirling engine was patented by Robert Stirling in 1816, and 
it was a scientifically remarkable power machine. Unfortunately, 
early Stirling engines were not competitive to IC engines mainly 
because of their small power-to-weight ratio and slow response 
to speed and power changes, and were literately ignored for quite 
a long time during the 19th and 20th centuries. Yet, their advanta¬ 
ges of high thermal efficiency, low maintenance requirement, flex¬ 
ibility on energy sources, and safe to operate have attracted 
significant interest in recent decades due to the world’s ever grow¬ 
ing energy demand, the depletion of fossil fuel, and the urgency to 
solve the worsening global warming problem. In particular, their 
ability to use solar energy has made them one of the most promis¬ 
ing solutions to global warming. Now Stirling engines have been 
the focus of many academic studies, both experimentally and 
numerically [1-3]. A new industry centered on Stirling engines 
and their applications has emerged [4]. Due to its high energy con¬ 
version efficiency, this 200-year-old invention even found its way 
to NASA’s deep space project [5], 

Stirling engines come in many different types; they are com¬ 
monly classified into three different catalogues according to their 
configurations, namely ot, p, and y types. Among these configura¬ 
tions, the p-type Stirling engines normally run at a high tempera¬ 
ture difference (HTD) between the hot and cold ends of the 
engine, making them very efficient. However, burning fossil fuel 
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or a complex solar energy collecting system is needed to sustain 
the high temperature at the hot end of the engine. The high tem¬ 
perature also gives rise to higher thermal radiation losses, higher 
thermal stress, and the requirement on expensive high-tempera- 
ture resisting materials. The y-type engines, on the other hand, 
normally run on a low or moderate temperature difference, en¬ 
abling them to be far more flexible on energy sources (biomass 
or waste heat) or requiring much simpler solar energy collecting 
techniques. In addition, the material requirement of such engines 
is less demanding due to the lower temperature at the hot end of 
the engine; and they can be built cheaply, which is ideal for domes¬ 
tic applications. However, the major disadvantage of the y-type en¬ 
gines is that they tend to be less efficient and less powerful than 
the p-type engines mainly because of low temperature difference. 

The studies on y-type engines are comparatively less than those 
on the P-type engines in open literature. Cinar and Karabult [6] de¬ 
signed and manufactured a small gamma-type Stirling engine with 
276 cc swept volume. The maximum power output of the engine 
was 128.3 W when charged with 4 bar helium and heated to 
1273 K at the hot end. Their work demonstrated that when run 
on HTD, a gamma-type Stirling engine is also capable of producing 
respectable power output. Parlak et al. [7] reported a thermody¬ 
namic analysis of a gamma-type Stirling engine using a quasi-stea¬ 
dy flow model which gave more precise results for the pressure 
profile than previous models. They concluded that a thermal 
efficiency of 25% is achievable for a gamma-type Stirling engine 
charged with nitrogen to 6.5 bar and a hot end temperature of 
873 K. Kongtragool and Wongwises [8] built and test a twin power 
piston and a four power piston y-type LTD Stirling engines 
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Nomenclature 



An 

wall area inside the expansion chamber 

W 

work (J) 

At2 

wall area inside the compression chamber 

W 

power (W) 

Cp 

constant pressure specific heat (kj kg -1 K _1 ) 

Y 

displacement (m) 

Cv 

constant volume specific heat (kj kg -1 K ') 



H 

convective heat transfer coefficient (W nr 2 K _1 ) 

Greeks 


h 

enthalpy per unit mass (kj kg 1 ) 

oc 

thermal diffusivity (m 2 s _1 ) 

k 

thermal conductivity (W nr 1 I<~ ] ) 

p 

rotation angle of displacer journal on the crankshaft 

U 

length of the power piston linkage bar (m) 

£ 

regeneration effectiveness 

h 

length of the power piston connection rod (m) 

p 

density (kg nr 3 ) 

h 

length of the displacer linkage bar (m) 

n 

engine efficiency 

U 

length of the displacer connection rod (m) 

e 

rotation angle of power piston journal on the crankshaft 

la 

height of displacer (m) 

o p 

phase angle between displacer and power piston 

m 

mass (kg) 

p 

fluid dynamic viscosity (kg nr 1 s _1 ) 

p 

pressure (N itT 2 ) 

CD 

engine speed (rad s~ ] ) 

Q 

heat transfer rate (W) 



R 

gas constant (kj kg -1 1< -1 ) 

Superscripts 

Rt 

outer radius of the power piston (m) 

n 

time-step number 

r 2 

inner radius of the displacer cylinder (m) 



Ra 

outer radius of displacer (m) 

Subscripts 

Rn 

thermal resistance of the expansion chamber (K W~') 

c 

compression chamber 

Rt 2 

thermal resistance of the compression chamber (I< W _1 ) 

d 

displacer 

n 

crank radius of the power piston (m) 

e 

expansion chamber 

r 2 

crank radius of the displacer (m) 

f 

end of a cycle 

T 

temperature (K) 

H 

hot end 

t 

time (s) 

i 

flow direction from regenerative channel to compres- 

t w 

thickness of displacer cylinder (m) 


sion chamber 

U 

internal energy (kj) 

j 

flow direction from regenerative channel to expansion 

U 

internal energy per unit mass (kj kg -1 ) 


chamber 

V 

volume (m 3 ) 

L 

cold end 

V 

velocity (m s -1 ) 

P 

power piston 


powered by simulated solar energy. The two engines’ performance 
were tested with air at atmospheric pressure; they yielded maxi¬ 
mum power and thermal efficiency of 11.8 W and 0.494%, respec¬ 
tively for the twin power piston engine and 32.7 W and 0.809%, 
respectively for the four power piston engine. Although both en¬ 
gines’ performance is not impressive, the author indicated that 
the engines can be improved by increasing the precision of engine 
parts and heat source efficiency and adopting helium or hydrogen 
as working fluid. The underwhelming performance of the engines 
could be further attributed to the fact that they were not optimized 
during the design stage. Karabulut et al. [9] pointed out that many 
parameters such as thermal and physical properties of working 
gas, temperature difference between hot and cold ends, charge 
pressure of working gas, and regenerator effectiveness affect the 
performance of a Stirling engine; and increasing or decreasing 
the values of these parameters often not only bring improvement 
on one aspect but also introduce negative impact on another. For 
example, increasing the heat transfer coefficient of the heater by 
adding more flow channels would also result in more flow losses 
and an increase on dead volume. These contradicting effects on en¬ 
gine performance highlight the need to perform parameter optimi¬ 
zation during the design stage. Cheng and Yu [10] argued that the 
design process of a new Stirling engine can be time consuming and 
expensive because the determination of geometrical parameters of 
the engine is merely based on the experience of the designers or 
insufficient experimental information. Therefore they developed 
a numerical model and applied it on a [3-type rhombic-drive Stir¬ 
ling engine. The model took non-isothermal effects, the effective¬ 
ness of the regenerative channel, and the thermal resistance of 
the heating head into account and was proved to be capable of 


identifying the most crucial geometrical parameters affecting the 
engine’s performance. However, many losses in real engines re¬ 
main unaccounted for in this model, hence its accuracy needs to 
be further verified by experimental data from real engines. 

Timoumi et al. [11] indicated that the actual Stirling engine effi¬ 
ciency remains very low compared to the high theoretical effi¬ 
ciency. Real Stirling engines involve very complex compressible 
fluid dynamics, thermodynamics, and heat transfer phenomena, 
and an accurate description and understanding of these highly 
transient phenomena and the engine losses is necessary to opti¬ 
mize engine’s design parameters. Several engine losses have been 
identified by Walker [12] as conduction losses in the heat exchang¬ 
ers, dissipation by pressure drop, shuttle and gas spring hysteresis 
losses, and according to Timoumi et al. [11], these losses are not 
usually accounted for in the published work due to their complex¬ 
ity. Costea et al. [13] also indicated that irreversibilities in the ther¬ 
modynamic cycle have significant importance in predicting the 
performance of Stirling engines, and analyses based on conven¬ 
tional entropy or exergy techniques do not relate the irreversibili¬ 
ties to the physical phenomena that cause them. Therefore, the 
objective of this study is to modify Cheng and Yu’s model to in¬ 
clude additional engine losses and then apply it on a charged twin 
power piston y-type Stirling engine to study some critical geomet¬ 
rical and operational parameters to the engine’s performance. 

2. Mathematical model 

We design and built a prototype air- or helium-charged, twin- 
power-piston y-type Stirling engine as illustrated in Fig. 1. In a 
charged Stirling engine, the entire engine, including the crank 
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system, is sealed inside a pressure vessel; therefore it is normally 
equipped with a built-in alternator to convert mechanical energy 
to electric energy, which can be readily drawn outside the pressure 
vessel via electric cables. If the engine is used to drive mechanical 
devices such as fans, pumps, and compressors, an electric motor is 
needed to convert electric energy back to mechanical energy. This 
practice introduces inevitable energy conversion losses. For exam¬ 
ple, even the efficiency of both alternator and electric motor 
reaches 80%; the two stages of energy conversion still result in al¬ 
most 40% of total energy loss. To avoid such substantial energy 
loss, a pair of magnetic couplings is employed here to transmit en¬ 
gine’s power through the wall of the pressure vessel to an external 
power shaft. The rotation of magnetic couplings would produce 
some eddy current on the metal casing; however, this energy loss 
is much lower than the energy conversion losses. This feature al¬ 
lows the direct use of the engine’s mechanical energy to drive 
mechanical devices. A schematic of the engine’s geometrical 
parameters is given in Fig. 2, where the engine has been divided 
into three control volumes: the expansion chamber, the regenera¬ 
tive channel, and the compression chamber. In this configuration, 
the walls of the expansion and compression chambers also serve 
as the heater and the cooler, respectively; therefore, there is no 
external heating or cooling channels. The heater and the cooler 
are respectively constructed with electric heating coils and fluid 
jacket surrounding the walls of expansion or compression cham¬ 
bers, thus the heating or cooling areas include not just the top or 
bottom plates but also parts of the side walls of the displacer cyl¬ 
inder. According to Fig. 2, the displacements of the piston and dis¬ 
placer, Yp(t) and Y d (t ) are: 

Yp(t) =-r, sin 0+yif^rfcos^e+i 2 (1) 

Y d (t) = —r 2 sin 0 + \Jl\ - r\ cos 2 p + U + k (2) 

where p = 0 - 0 P , 0 P is the phase angle between the piston and dis¬ 
placer. With the determination of V p (t) and Y d (t), the volumes of the 
expansion and compression chambers can be expressed as: 



Fig. 2. Definition of engine geometrical parameters. 


V e (t) = (/ci - Y d )nR 2 2 

(3) 

V c (t) = ( Y d -l d - l c2 )nRl + 2(la - Y p )nR 2 

(4) 
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where the factor of 2 multiplied to the second term of Eq. (4) ac¬ 
counts for the twin power piston configuration. The velocities of 
piston and displacer are then calculated by respectively differentiat- 


ing Eqs. (1) and (2): 


. . „ r?m cosd sin 9 

v p (t) = -r-ita cos 9+ ' 

\Jl 2 - t \ cos 2 8 

(5) 

, , „ rococos/I sin/! 

Vd{t) = -r 2 to cos p + 2 .= 

V*3 - r 2 COS 2 P 

(6) 

By using ideal-gas equation of state, the pressure in the expan¬ 
sion and compression chambers are: 

D m e RT e D m c RT c 

‘ e ~ V e ’ ‘ L V c 

(7) 

Then the energy conservation law is applied respectively to the 
expansion and compression chambers and regenerative channel; 
this read: 

dll . . / v 2 \ 

-gp = Q-in.e - Qloss ~ W mt<e - ifl + yj 

(8) 

— Qin,r — tTl(hj: n — h out ) r — — [y\ n — V o Ut ) r 

(9) 

dll ■ ■ ■ ( v 2 \ 

= Qin,c + Qioss — W ou t,c + rn yh + J 

(10) 


where Q j0SS represents the heat losses from the expansion chamber 
to the compression chamber. According to Timoumi et al. [11], the 
dominant heat losses are internal conduction heat loss and shuttle 
heat loss; hence the term is approximated by: 


Q/oss — Q.cond “E Qshtl (11) 

The internal conduction heat loss is the energy loss due to inter¬ 
nal conduction between the hot and cold parts of the engine. In the 
current configuration, this loss reads: 


Qcond — 


27iR 2 t w (TH - Ti) 
Lnnd 


( 12 ) 


{ h = h e , v = v e \ expansion 
h in = h e ,h oul: = hi, Vin = v e , v 0M = vc, regenerative 
h = h,, v = vc compression 

(14) 


{ h = hj, v = Vj \expansion 
h in = h c ,h out = hj, v in = V c , Vout = Vj \regenerative 
h = h c , v = v c ', compression 

(15) 


where subscripts i and j represent the flow properties from the 
regenerative channel to the compression chamber and the flow 
from the regenerative channel to the expansion chamber, respec¬ 
tively. For flow passing the regenerative channel, the temperatures 
Tj and 7) are approximated by: 

Ti = T e + e(T c - T e ), m > 0 

(16) 

Tj = Tc + e(T e -Tc), m <0 


where e is the effectiveness of the regenerative channel which poses 
decisive impact on the performance of the Stirling engine. The 
internal energy and enthalpy of the working fluid can be calculated 
by: 

u = c v T, h = c p T (17) 


The heat transfer rates in the expansion and compression cham¬ 
bers are approximated by: 



(18) 


where T H , R n and T L , R t2 are the temperatures and thermal resis¬ 
tances respectively at the walls of the expansion and compression 
chambers. In Cheng and Yu [10], the values of R n and R t2 are fixed. 
In reality, as the heat transfer areas of the expansion and compres¬ 
sion chambers vary with time, these values should also vary with 
time. The heat transfer areas of the expansion and compression 
chambers A n and A a can be written as: 


An = nR 2 + 2nR 2 (l c t - Y d ), 

Aa = (nR 2 2 - 2nR]) + 2nR 2 (Y d - l c2 - l d ) (19) 


where t and l cond are the thickness of the displacer cylinder and 
the distance between the hot and cold ends of the engine, 
respectively. The shuttle heat loss is rooted in the reciprocating 
motion of the displacer, which absorbs a quantity of heat from 
the hot end and restores it at the cold end. Chang and Baik 
[14] derived a simple expression to account for the shuttle heat 
loss as: 



where af is the thermal diffusivity of the working gas, and k s the 
thermal conductivity of solid material. The above equation is sim¬ 
plified from Chang and Baik’s original equation by assuming the 
displacer cylinder wall and the displacer are constructed by the 
same solid material. In Eqs. (8)-(10), we assume that the mass 
flow out of the expansion chamber is positive. Since the working 
fluid flows in and out of the expansion and compression cham¬ 
bers repeatedly via passing the regenerative channel, a convection 
scheme is needed to determine the fluid properties flowing in 
and out of each control volume. Here, a simple upwind scheme 
is adopted; that is: 


and R n and R t2 become: 


Rn 


1 

TWt' 


Ra = 


HA t: 


( 20 ) 


where H is the convective heat transfer coefficient on the walls of 
expansion and compression chambers. Karabulut et al. [15] sug¬ 
gested the values of H for air and helium in their engine are 
447 W itt 2 IC 1 and 2392 W m~ 2 K -1 , respectively. The value for 
helium will be tested in the current study. The power output is cal¬ 
culated by: 


■ Pe(K +1 -0 Pc(V n c +1 -V n c ) 

out,e — - 7T- ■ vv out,c — - 


At 


At 


( 21 ) 


where superscript n represents the marching time step number. 
Substituting Eqs. (12)-(15) into Eqs. (8) and (10), the temperatures 
in the expansion and compression chambers can be obtained: 


r"+i 


= 1 - 


At 


rh„c„R, 




At 


m e c v 


e L y Afi 

Th 

Rn “ 


yfl 


PeiV^ - V") 

At 


- m 




- Qto 


( 22 ) 


for m > 0, or 
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r „+, = M 


At 


m e c v R t i 


I’ll 


At 

-F -=- 

m e c v 

for m < 0; and 

rn 1 


Th Pe(V n e +1 - V n e ) 
Rn 


At 


- m ( hj + -j -) - Q (0 , 


r" + = 1 - 


At 


m c c v R a 


■j’H 


+ — 


At 
m c c v 

for rh > 0, or 
T" +1 = ( 1 - 


r L p c (v c n+1 - v c ") 

Rt 2 


At 


+ rn ( hj + ) + Qjoss 


At 


m c c„R t2 




+ — 


At 


HtrC,, 


Ti. p c (v n c +1 - v c n ) 

Rn At 


+ rn ( h c + ^ ) + Q| 0 , 


(23) 


(24) 


The magnitude of frictional losses depends on many factors in man¬ 
ufacturing the engine such as the accuracy of machining, clearance 
between the piston and cylinder, and the alignment of engine parts, 
especially the crank assembly. Misalignment of crank assembly 
would inadequately compress shaft bearings and cause a dramatic 
rise of bearing friction. The friction loss magnitude cannot be pre¬ 
dicted; however, it is reasonable to assume that friction loss is con¬ 
stant in every engine cycle, and the power lost to friction Wf can be 
approximated by: 


W f = ^W fc (30) 

Finally, the thermal efficiency of the engine can be obtained by: 


W„, 


Jo l Q.m,dt 


(31) 


(25) 


3. Results and discussion 


for rh < 0, where the bar over a quantity represents the averaged 
value of that quantity between time steps n and n + 1: 




(26) 


In Eqs. (22)-(25), the wall temperatures of the expansion and com¬ 
pression chambers are assumed to be T H and T L , respectively. In 
reality, however, T H or T L is maintained at the external walls of 
the expansion or compression chambers, and heat is transferred 
into or out of the expansion or compression chambers from the in¬ 
ner walls of both chambers via mechanism of conjugate heat trans¬ 
fer. In this case, the inner wall temperature of the expansion 
chamber is lower than T H , and that of the compression chamber is 
higher than T L . To take the mechanism of conjugate heat transfer 
into account, the wall temperatures of the expansion and compres¬ 
sion chambers T H and T L should be respectively modified as: 


Twh = 


;T H + HT e 

t + H 

LW 


Twl = ' 


'-t l + ht c 


H 


(27) 


In numerical model, the regenerator is assumed at the annular 
between the displacer piston and cylinder. The mass flow rate rii 
is calculated by solving the viscous velocity profile in the annular 
regenerative channel and integrating the mass flow rate across 
the channel; this yields: 


m = - 


pn 

8p 




In (R 2 /R d ) 


2npvd 

~ In (R 2 /Rd) 
npv d 

In(R 2 /Rd) 


^ (r] In R 2 - R 2 d ln R d J 
In R 2 (R 2 2 -R 2 d ) 



(28) 


In the real engine, the regenerator includes not only the annular 
but also the space inside the displacer piston where very thin cop¬ 
per wires have been put inside to serve as the regenerative mate¬ 
rial. This is done by drilling a number of small holes on the top and 
bottom plates of the displacer piston, allowing working gas to pass 
through the interior of the displacer piston. Hence, the approxima¬ 
tion of regeneration channel adopted in the model might be a 
source of error. The network output per cycle is calculated by the 
following integration: 

W out = f’ P e dV e + [' P c dV c ~ W f (29) 

Jo Jo 


where t/ represents the time at the end of one cycle, and Wf ac¬ 
counts for frictional losses from all moving parts of the engine. 


The values of the geometrical parameters of the baseline Stir¬ 
ling engine are given in Table 1. The whole engine, including the 
crank system and fly wheels, is housed inside a well-sealed 3- 
rnm thick stainless steel pressure vessel, allowing the engine to 
be charged up to 10 bars with air or helium. Since the current en¬ 
gine has two power pistons and one displacer piston, the crank 
shaft is consisted of six rectangular sections and three connection 
rod pins, all joined by welding to provide superior strength. Ex¬ 
treme care was taken to minimize misalignment of these pieces 
to reduce inadequate compression on shaft bearings. The body of 
the power piston was made of aluminum, and that of the power 
cylinder was of brass. The clearance between the power piston 
and cylinder is 0.02 mm, and the contact surfaces of both were pol¬ 
ished to super-finish quality. The displacer piston was constructed 
using two 2-mm thick circular aluminum top and bottom plates 
covered with 0.1-mm stainless steel sheet to form the lateral skin. 
The displacer cylinder was made of steel pipe with inner diameter 
of 160 mm and wall thickness of 3 mm. The cooler is a water jacket 
welded around the upper part of the displacer cylinder, through 
which fresh cooling water is circulating to maintain low tempera¬ 
ture, and the heater is a thermostat controlled electric heater, 
capable of heating up to 873 K, wrapped on the lower part of the 
displacer cylinder. Fig. 3a shows the entire engine, where the pres¬ 
sure vessel and the external magnet of the magnet couplings can 
be seen. In Fig. 3b, the interior of the pressure vessel is shown, 
including the crank system, twin power cylinders, and the fly¬ 
wheels. A torque transducer, KISTLER 4520A100 with accuracy of 
0.1 N-m, was installed between the engine’s power shaft and the 
loading to measure the engine’s torque. A photo tachometer, TEC- 
PEL UT371 with accuracy of 0.01 RPM, was employed to measure 


Table 1 

The values of the geometrical parameters of the baseline Stirling engine. 


R, (m) 
R 2 (m) 
Rd(m) 
r, (m) 
r 2 (m) 
h (m) 

1 2 (m) 

1 3 (m) 
U (m) 
Id (m) 
la (m) 
l C 2 (m) 
t(m) 
Bp (°) 


0.0256 

0.079 

0.077 

0.025 

0.025 

0.13 

0.045 

0.095 

0.155 

0.148 

0.425 

0.221 

0.003 
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(b) 

Fig. 3. The prototype Stirling engine, (a) the pressure vessel, and (b) the interior of 
the pressure vessel, including twin power cylinders, crank system, and flywheels. 


the engine’s rotation speed. The measurements from the torque 
transducer and the tachometer allow the engine’s power to be cal¬ 
culated. Since using a real Stirling engine to conduct a parametric 
study is not practical, in this study, the prototype Stirling engine 
only serves the purpose of providing data to calibrate the numeri¬ 
cal model. 

The numerical code is first validated using the results in Cheng 
and Yu [10], where a (3-type rhombic-drive Stirling engine was 
simulated, and results in terms of volume, temperature, pressure, 
and heat transfer variations in Stirling cycles were well docu¬ 
mented for a base-line engine. Details of the geometric and opera¬ 
tion parameters of the base-line case can be found in [10] and will 
not be repeated here. The calculation in [10] did not take heat- 
transfer and friction losses, and effect of conjugate heat transfer 
into account; therefore, these losses are set to zero, and conjugate 
heat transfer effect is ignored in the current code for code-valida¬ 
tion purpose. As shown in Fig. 4, the comparison of p-v diagram of 
the expansion and compression chambers indicates good agree¬ 
ment between the results obtained by the current code and those 
in [10]. The calculated engine power of 16.7 W is also very close to 
the value of 16.75 W reported in [10]. This exercise has verified the 
correctness of the current numerical code. 


Expansion 



Fig. 4. Comparison of expansion and compression chambers p-v diagrams of the 
current code and Cheng and Yu [10]. 


Before the numerical simulation of the gamma-Stirling engine, 
the value of Wf c in Eq. (29) needs to be determined. Since this value 
is highly dependent on the manufacture quality of an engine, it can 
only be determined by measurement. A simple method is used to 
estimate this quantity. First, the prototype engine was run with 
air under atmospheric pressure, and T H and T L were fixed at 
473 K and 300 K, respectively. In these conditions, the engine only 
yielded small power. A small loading was applied, and the mea¬ 
sured engine power was 3.6 W at the engine speed of 104 rpm. 
Then the loading was removed, and the engine speed was in¬ 
creased by 196 rpm to 300 rpm. When the loading was removed, 
the engine’s power was consumed entirely by friction loss. Assum¬ 
ing that engine’s power output is the same at both engine speeds, 
the power output of 3.6 W is equivalent to W f , and the increased 
engine speed of 196 rpm is equivalent to the rotational speed. In 
this case, the value of Wf c is 1.08J/cycle. This method only gives 
a rough estimation on Wf c because the engine’s power output actu¬ 
ally depends on engine’s speed. 

The thermodynamic cycle of the prototype engine under a set of 
baseline operation conditions is simulated by the code. The base¬ 
line operation parameters are as follows: engine speed =-400 
rpm, hot end temperature T H = 700 K, cold end temperature 
T l = 300 K, the working gas is helium charged to 2 bar, convective 
heat transfer coefficient at hot and cold ends H = 2396 W m~ 2 K ’, 
and regeneration effectiveness e = 0.3. Here, the engine speed is 
typical for a y-type Stirling engine, and the predefined hot end 
temperature is reachable by burning biomass or using simple 
one-axis tracking parabolic trough solar collector. Figs. 5-7 respec¬ 
tively show the volume variations, temperature variations, and p-v 
diagrams of the expansion and compression chambers. Fig. 5 
shows that the variations of volume in the expansion and compres¬ 
sion chambers are similar; however, the minimum compression 
volume is larger than the minimum expansion volume because 
there are dead volumes in the twin power cylinders. As can be seen 
in Fig. 6, the variation of temperature in the expansion chamber is 
about 70 K, and that of the compression chamber is about 50 K. 
Cheng and Yu [10] reported the simulated maximum expansion 
chamber temperature is a bit higher than the hot end temperature, 
and the minimum compression chamber temperature is a bit lower 
than the cold end temperature. They attributed the phenomena to 
the high compression ratio of their [3-type engine. Such overshoot¬ 
ing phenomena are not observed here presumably due to the much 
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Fig. 5. Volume variations in the expansion and compression chambers under 
baseline conditions. 



Fig. 6. Temperature variations in the expansion and compression chambers under 
baseline conditions. 


lower compression ratio in the current engine and the introduction 
of the conjugate heat transfer effect in the numerical model. From 
Fig. 7, it can be seen that during one engine cycle, the expansion 
chamber yields positive net work, whereas the compression cham¬ 
ber produces negative net work. The absolute value of the positive 
net work is larger than that of the negative net work, giving rise to 
positive work output to the engine’s shaft. Therefore, engine’s effi¬ 
ciency can be improved by promoting the positive net work in the 
expansion chamber and reducing the negative net work in the 
compression chamber. The calculated engine power is 159.77 W 
for the baseline conditions. This value, however, is significantly 
overestimated when compared with the measured power of 
102 W. The calculated heat transfer rate into the expansion cham¬ 
ber is 2.3 kW, which is close to the input electric power. Some of 
the reasons contributing to the overestimation of engine power 
could be that the values of convective heat transfer coefficient H 
and the regeneration effectiveness e are too high in the model. 
By resetting the values of H and e to 1500 Wm^K -1 and 0.2, 
respectively, the calculated engine power becomes 146.1 W, much 


Expansion 



Fig. 7. p-v diagrams in the expansion and compression chambers under baseline 
conditions. 

closer to the measured value. Other reasons responsible for the 
inaccuracy of the model could be more fundamental, such as the 
assumption of annular regeneration channel and uniform pressure 
inside the expansion and compression chamber at any given mo¬ 
ment. To take the effects of these factors into account, a 3D CFD 
simulation using moving grid technique and solving Navier-Stokes 
equations is needed. Yet, such undertaking is very sophisticated 
and time consuming and is not appropriate for parametric study. 
In the following calculations, the smaller H value will be adopted. 

Next, the effects of four important geometrical and operational 
parameters, namely the regeneration effectiveness, the crank ra¬ 
dius of the power piston the initial pressure of working gas, 
and the rotation speed, on the engine performance are examined. 
In this analysis, the working gas is helium, and the hot- and low- 
end temperatures are fixed at 700 K and 300 K, respectively. The 
regeneration effectiveness has been well known as one of the most 
influential factors promoting the engine’s power and efficiency. 
The value of r, determines the swept volume of the power pistons, 
which, according to [6], poses crucial impact on the engine’s 



Fig. 8. The effects of regeneration effectiveness on engine's power and efficiency. 
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Fig. 9. The effects of power piston crank radius on engine's power and efficiency. 

performance. The initial pressure dictates the amount of working 
gas inside the engine and is also influential to power output. The 
engine speed is yet another very crucial factor on its power and 
efficiency. The effects of these parameters are respectively shown 
in Figs. 8-11. In Fig. 8, the values of r t , initial pressure, and engine 
speed are fixed at 0.025 m, 2 bars, and 400 rpm respectively. The 
results clearly indicate that both the engine’s power and efficiency 
are proportional to e. Its effect is more pronounced on the engine’s 
efficiency which increases more than three times when e increases 
from 0.05 to 0.5. To investigate the effect of the values of e, ini¬ 
tial pressure, and engine speed are fixed at 0.3, 2 bars, and 
400 rpm, respectively. It is noticeable in Fig. 9 that both engine’s 
power and efficiency increase almost linearly with the increase 
of ri; especially the power is increased by a factor of 2 as ri changes 
from 0.02 m to 0.04 m. In terms of investigating the effect of initial 
pressure, the values of e, it, and engine speed are fixed at 0.3, 
0.025 m and 400 rpm, respectively. As can be seen in Fig. 10, the 
power and efficiency also increases almost linearly as pressure in¬ 
creases. This parameter is very effective on both power and effi¬ 
ciency. Finally, to examine the effects of engine speed, the values 
of e, ri, and initial pressure are fixed at 0.3, 0.025 m, and 2 bar, 



Fig. 10. The effects of initial change pressure on engine’s power and efficiency. 



Fig. 11. The effects of engine speed on engine’s power and efficiency. 

respectively. Fig. 11 indicates that engine speed is of great impor¬ 
tance on engine power. As the engine speed increases from 
100 rpm to 400 rpm, the power increases from 40 W to 160 W, a 
leap of four folds. This seems to suggest that engine power can 
be readily boosted by just increasing engine speed. However, the 
y-type engine tends to have a large displacer piston. Particular in 
the current configuration, the displacer also houses the regenera¬ 
tion material, making it a heavy component. Its mass is between 
0.5 and 0.8 kg, depending on the amount of regeneration material 
put inside. The heavy weight of the displacer limits the engine 
speed because its reciprocation motion inside the displacer cylin¬ 
der results in huge change of momentum at high speed. The 
momentum change acting on the crank shaft causes observable 
unbalance and vibration as engine speed reaches 400 rpm. There¬ 
fore, this type of engine is not engineered to run at high speed, 
and the engine speed of 400 rpm is more or less the safe top speed 
of the current y-type Stirling engine. The impact of engine speed 
on efficiency is less dramatic as efficiency only increases from 
0.55 to 0.79 when engine speed increases from 100 to 400 rpm. 

4. Conclusions 

In this paper, a parametric study has been conducted on a he¬ 
lium-changed twin-power-piston y-type Stirling engine. A baseline 
engine has been constructed to correct some parameters in the 
numerical model. Four important geometrical and operational 
parameters have been investigated, including the regeneration 
effectiveness, the crank radius of the power piston, the initial pres¬ 
sure of working gas, and the rotation speed. It has been found that 
regeneration effectiveness poses the most prominent effect on effi¬ 
ciency, while increasing engine speed is most effective on the en¬ 
gine power within the range of engine speed investigated. The 
crank radius of power piston and initial gas pressure, on the other 
hand, are also very influential on engine power. Their effects on en¬ 
gine efficiency, however, are not as diametric. This study offers 
invaluable information on the optimization of the prototype y-type 
Stirling engine and the design of other y-type engines with similar 
construction. 
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